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ABSTRACT 

Context. Debris disks are commonly considered to be a by-product of planet formation. Structures in debris disks induced by planet- 
disk interaction are promising to provide valuable constraints on the existence and properties of embedded planets. 
Aims. We investigate the observability of structures in debris disks induced by planet-disk interaction with future facilities in a 
systematic way. High- sensitivity, high angular resolution observations with large (sub-)mm interferometers and large space-based 
telescopes operating in the near- to mid-infrared wavelength range are considered. 

Methods. The observability of debris disks with the Atacama Large Millimeter/submillimeter Array (ALMA) is studied on the basis of 
a simple analytical disk model. Furthermore, N-bo&y simulations are used to model the spatial dust distribution in debris disks under 
the influence of planet-disk interaction. From these simulations, images at optical scattered light to millimeter thermal re-emission are 
computed. Available information about the expected capabilities of ALMA and the James Webb Space Telescope (JWST) are used to 
investigate the observability of characteristic disk structures with these facilities through spatially resolved imaging. 
Results. Our simulations show that planet-disk interaction can result in prominent structures in the whole considered wavelength 
range. The exact result depends on the configuration of the planet-disk system and on the observing wavelength which provides 
the opportunity of detecting and characterizing extrasolar planets in a range of masses and radial distances from the star that is not 
accessible to other techniques. Facilities that will be available in the near future at both considered wavelength ranges are shown to 
provide the capabilities to spatially resolve and characterize structures in debris disks that arise because of planet-disk interaction. 
Limitations are revealed and suggestions for possible instrument setups and observing strategies are given. In particular, ALMA is 
limited by its sensitivity to surface brightness, which requires a trade-off" between sensitivity and spatial resolution. Space-based mid- 
infrared observations will be able to detect and spatially resolve regions in debris disks even at a distance of several tens of AU from 
the star, where the emission from debris disks in this wavelength range is expected to be low. 

Conclusions. Both ALMA and the planned space-based near- to mid-infrared telescopes will provide unprecedented capabilities to 
study planet-disk interaction in debris disks. In particular, a combination of observations at both wavelengths will provide very strong 
constraints on the planetary/planetesimal systems. 

Key words. Techniques: high angular resolution - Techniques: interferometric - Planets and satellites: detection - Planet-disk inter- 
action - Infrared: planetary systems - Submillimeter: planetary systems 



1. Introduction 

The dust detected in debris disks is thought to be removed from 
those systems by the stellar radiation on time scales that are short 
compared to their ages. This means that the dust must be tran- 
sient, or more likely continuously replenished by ongoing colli- 
sions of bigger objects such as planetesimals left over from the 
plane t formation process (for a recent review see, e.g., | Krivov| 
|2010| ). The presence of planets and debris disks is hence thought 
to be correlated. In a system with a debris disk and one or more 
planets, one would expect gravitational interaction between the 



dust grains and the planet, trapping them into resonan ce ( |Wyatt 
[20061 1 Wolf et al.|2Q07l [Stark & Kuchner|2008| [2009) ). This re- 
sults in structures in the disk that may be observable, which in 
turn provides a method to infer and characterize planets that are 
in a regime of masses, brightnesses, and radial distances from the 
star that is not accessible through other techniqu es such as radial 
velocity measure ments or direct imaging ( e.g. 
[20081 [MaroiseFa l. 2008; Kalas et al.|2008| ). 

Clumpy structures in debris disks have b een obser ved in sev- 
eral cases (e.g., 6 Eri, [Greaves et al.[p^98l AU Mic, |Liu||2004| 



Udry & Mayor 



HD 107146, [Corder et al.[[2QQ9l [Hughes et al1[2QTT] ). In this 
work, we investigate the observability of structures in debris 
disks in a systematic way. We set up several initial conditions for 
the planetary mass and orbit as well as for the dust distribution in 
the disk from which we simulate the spatial dust distribution us- 
ing A/^-body simulations. Our results from th e dynamical simula- 
tions are consistent with earlier work s (e.g., Holmes et al.|2003( 
|Wyatt|2006"l [Stark & Kuchner|20 08'). In contrast to these works 
which focused mostly on the phenomenological and theoretical 
description of the structures as well as on the processes of planet- 
disk interaction, the goals of the present work are the evaluation 
of the observability and the development of strategies for obser- 
vations of planet disk interaction. Thus, we use the A/^-body sim- 
ulations as a tool to produce realistic structures in debris disks, 
but discuss these structures only briefly emphasizing on the con- 
sequences of difl'erent planet-disk configurations on the observ- 
ability. To make predictions on the feasibility to detect and spa- 
tially resolve characteristic structures, available information on 
the capabilities of the Atacama Large Millimeter/submillimeter 
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Array (ALMA) and James Webb Space Telescope (JWST) are 
used representative for near- future facilities. 

We describe the approach used in our N-ho&y simulations 
and for the image generation in Sect. [2] The results from the 
dynamical simulations are presented and briefly discussed in 
Sect. [3] In Sects. |4] to [6j we evaluate how large-scale disk struc- 
tures may be observed with difl'erent facilities. Conclusions are 
drawn in SectH 



power-law with exp onent -3.5 as exp ected from an equilibrium 
collisional cascade ( |DohnanyTl|1969 ). This approach is used to 
mimic dust production through collisions of the parent bodies. 
Note that the local grain size distribution may be significantly 
diff'erent due to the redistribution of the grains through dynam- 
ical interaction with the planet as well as throug h the eff'ect of 
radiation pressure and Poynting-Robertson drag ( Moro-Martin 
,& Malhotra,2QQ2,.2QQ3j . 



2. Modeling planet-disk interaction in debris disks 2.1 .2. Grain size distribution 



In this section, we describe the initial conditions employed for 
modeling planet-disk interaction. We simulate the dynamical 
evolution of a large ensemble of dust particles in the gravita- 
tional potential of a central star under the influences of planetary 
perturbations, radiation pressure, Poynting-Robertson efl'ect, 
and stellar-wind drag. For these simulations, the code MODUST 
(|Rodmannl2006 ) is employed. Electromagnetic (Lorentz) forces 
on c harged dust particles are neglected ( Gustafson 1994; H olmes 
|et al.||2003| ). Mutual collisions of dust particles are not consid- 
ered in the dynamical modeling. A brief description of the meth- 
ods and numerical approach to simulate the dynamical evolution 
of the dust grains as well as of the approach to simulate images 
from the resulting density distributions can be found in the ap- 
pendix. 

2.1. Initial conditions 

We use several initial conditions (planetary orbit and mass as 
well as initial dust distribution) selected by analogy to our solar 
system and to other known debris disks in order to study the 
structures that can be expected depending on the configuration of 
the respective system. Since the simulation of the resulting dust 
distribution is very time-consuming, we are not able to perform 
a dense sampling of a large, high dimensional parameter space. 
Therefore, we focus on the following parameter regime: 

- One central star (no multiple systems), 

- a solar-type star as central star, photospheric emission re- 
alized by a black body with Teff = 5778 K, = 1.0 L©, 

= l.OMo, ^ = 0.35 ( |Gustafson|1994l ), 

- one planet (considered to dominate the dynamics of the sys- 
tem), 

- one initial dust disk (no initial multi-ring systems), 

- planetary orbit and initial dust disk are coplanar. 

Astronomical silicate with a bulk density o f 2.7g/cm^ 
( Draine & Lee|1984{ Weingartner & Draine| 2001) is employed 
for the chemical composition of the dust. In the following, the 
initial conditions for the dust distribution (spatial distribution 
and grain size distribution) and for the planetary orbit are de- 
scribed and motivated. 



2.1 .1 . Treatment of dust creation 

We assume the dust to be produced through collisions in a disk 
of planetesimals. It is then redistributed because of the eff'ects of 
stellar radiation and wind, and the gravitational interaction with 
the planet (Eq. |A.l| ). This is realized by placing the initial dust 
distribution at the same position as the planetesimal disk postu- 
lated to produce the dust. Particles that are lost because they are 
sublimated or ejected from the system are replaced by new ones 
from this reservoir of initial dust grains. Furthermore, the grain 
size distribution of all grains in the system applied follows a 



The grain size can be expressed by the corresponding value of /3. 
The correlation between these two quantities is given by 
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where a is the grain size, Q^r is the efficiency of radiation pres- 
sure on a grain, and are the luminosity and mass of the 
star, c is the speed of light, G is the gravitational constant, and 
p is the bulk density of the dust grains which are assumed to be 
spherical an d compact. The approximation was given by |Wyatt| 



et al. (1999]) and is valid in the geometrical optics approxima- 



tion. Since the critical parameter for the computations is j3, it is 
used in MODUST directly to describe the grain properties (instead 
of grain size and optical properties). Around a solar-type star, the 
approximation gives good results for jS < 0.5 ( Rodmann^2006) . 
The lower boundary of the applied grain size distribution is de- 
fined hy JS = 0.5 (a ^ 0.43 yum), because smaller grains (higher 
values of fi) are ejected from the system by radiation pressure. 
The upper boundary of the distribution is set to a grain size of 
2 mm (J3 ^ 1.06 X 10~^ for the applied chemical composition). 
The contribution of larger grains to the simulated observations is 
neglected. 

The grain size distribution is sampled with 50 grain sizes 
between the lower and the upper grain size, distributed logarith- 
mically (to properly sample the thermal re-emission of diff'erent 
grain sizes for the final images). For each size bin, the dust spa- 
tial distribution is simulated. A total number of 1000 test parti- 
cles is used in each size bin. The system evolves over five times 
the corresponding Poynting-Robertson time scale of grains of 
the corresponding ft at the outer edge of the disk, but not longer 
than the assu med ag e of the system (see motivation of the sin- 
gle runs. Sect. 2.1.4 ). At the end of each run, 100 snapshots of 
the particle distribution (Ro dmann^2006 ) are taken in the refer- 
ence frame corotating with the planet to virtually increase the 
number of particles in each size bin. The snapshots are equally 
distributed over the last 10% of each run time (see description of 
the diff'erent runs, Sect. |2.1.4| ). 

2.1 .3. Spatial dust distribution 

Gaussian distributions are applied to set the initial eccentricities 
and inclinations of the test particles. At the start of the integra- 
tion, the majority of dust grains have orbital eccentricities be- 
tween and ~ 0.4 (rms of the Gaussian distribution = 0.15). 
The full width at half maximum of the inclination distribution 
is set to 10°. These values are consistent with measurements on 



dust parent bodies in the Kuipe r Belt of the solar system ( [Jewitt 
et al.|1996 Vitense et al.|2010 ). The semi-major axes of the par- 
ticles are distributed following a power-law ( Wolf & Hillenbrand 
2003|) from a lowest to a highest value considered as the inner 



and outer radius of the disk. 
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2.1 .4. Explored parameter space 



Table 1. Initial conditions for the different runs with MODUST 



The parameter space explored is motivated by analogy to our so- 
lar system and to other known debris disks. We use the position 
and shape of the dust disk as well as the position, mass, and ec- 
centricity of the planet as parameters for our initial conditions. 
The different combinations of values explored for these parame- 
ters are listed in Table [T] and described and motivated below. 

A first sequence of initial conditions (la to Id) places the planet 
in 1:1 or in 2:1 resonance with the initial planetesimal belt. The 
initial dust disk is a narrow ring (7?out = 1-1 ^in) with constant 
surface density and two realizations of Ri^ (5AU for runs la 
and lb and 50 AU for runs Ic and Id). The planet has a mass of 
1 Mj and is on a circular orbit. The two assumed positions of the 
dust ring are motivated by our solar system (Asteroid Belt and 
Kuiper Belt) and typical positions of the dust i n other known 
debris disks around Sun-like stars (e.g., 6Eri, Backman et al. 
[2009; HD 105). An age (run time of the simulation) of 50Myr 
is assumed for th ese runs (e.g., the approximate age of HD 105; 
|Apai et al.|2QQ8l ). 

The second sequence of initial conditions (Ila to Ild) is moti- 
vated byjhe HD 107146 debris disk ( [Ardila et al.]|2nQ4l |ErteT| 
et al. 201 1| ). This disk is particularly well-suited for studying 



structures in debris disks because of its nearly face-on orien- 
tation and the large amount of complementary data available. 
Furthermore, (sub-)mm observations of the disk have been pub- 
lished by [Corder eFaT] ( [20091 ) and [Hughes et al.[ ( [20n] ), show- 
ing clumpy structures in the surface brightness that could be the 
signposts of planet-disk interactioi£] Here, we employ a set of 
initial conditions where the planet is orbiting at the inner edge 
of an extended disk with a large inner hole. The inner and outer 
radius of the disk are set to Ri^ = 70 AU and /?out = 250 AU. The 
radial surface density distribution is described by a power-law 
with an index of y = -0.5 (e.g., [Wolf & Hillenbrand 2003). This 
model is similar to the results of deta iled multi-wavele ngth mod- 
eling of the HD 107146 debris disk ( [Ertel et al.|2QTT] ), although 
a more complex radial density distribution has been found there. 
However, such a complex distribution might be the result of the 
interaction of the disk with a possible planet and the temporal 
evolution of the disk. We explore different parameters for the 
planet. It is on a circular orbit and has a mass of 0.5 Mj, 1.0 Mj, 
and 5.0 Mj in the runs Ila, lib, and lie. It is on an eccentric or- 
bit (e = 0.1) and has a mass of 1.0 Mj in run lid. The age of 
the system is increased compared to that employed in our se- 
quence I to account for the slightly higher age of HD 107146 
(e.g., [Roccatagliata et al.|2009| ). An age of 100 Myr is employed. 

A third sequence of models (Ilia to Illd) places the planet within 
a broad disk. Therefore we modify the initial conditions of our 
sequence II, so that the disk is now closer to the star, ranging 
from 7?in = 35 AU to 7?out = 210 AU. The inner radius is com- 
parable to that of the Kuiper Belt. The power-law index of the 
radial dust distribution and the parameters of the planet remain 
the same. 



2.2. Image creation 

We employ Mie theory to com pute the optic al properties of the 
dust using the tool mi ex ( Wolf & Voshchinnikov 2004). To ere- 



Run 


D r A T n 
Kin LAUJ 


D r A T Tl 

Rout LAUJ 


a 


Mpi LMjJ 


^ r A T Tl 

api LAUJ 




la 


5.0 


5.5 


0.0 


1.0 


5.0 


0.0 


lb 


5.0 


5.5 


0.0 


1.0 


3.15 


0.0 


Ic 


50.0 


55.0 


0.0 


1.0 


50.0 


0.0 


la 


jU.U 


jj.U 


U.U 


i.U 


J L.J 


U.U 


Ila 


70.0 


250.0 


-0.5 


0.5 


70.0 


0.0 


lib 


70.0 


250.0 


-0.5 


1.0 


70.0 


0.0 


lie 


70.0 


250.0 


-0.5 


5.0 


70.0 


0.0 


lid 


70.0 


250.0 


-0.5 


1.0 


70.0 


0.1 


llla 


35.0 


210.0 


-0.5 


0.5 


70.0 


0.0 


Illb 


35.0 


210.0 


-0.5 


1.0 


70.0 


0.0 


Illc 


35.0 


210.0 


-0.5 


5.0 


70.0 


0.0 


Illd 


35.0 


210.0 


-0.5 


1.0 


70.0 


0.1 



^ However, it is doubtful whether these structures are real, since they 
are only visible at a signal-to-noise ratio (S/N) of ^ 5 to 6 (2 to 3 above 
the average surf ace brightness of the dis k) and their position i n the two 
observations by |Corder et al.| ( [2009] ) and |Hughes et al.|f2011| ) seems to 
be inconsistent. 



ate images from the derived dust distributions, we compute the 
thermal re-emission and the scattered stellar light for each dust 
particle at a set of observing wavelengths. The fluxes are then 
projected onto the plane of the sky considering the scattering 
phases. This can be performed for arbitrary orientations of the 
disk (assuming an optically thin debris disk). The image exten- 
sion in one dimension used is 301 pixels for sequence I (resulting 
in a pixel resolution of 0.1 AU for runs la and lb and 1.0 AU for 
runs Ic and Id) and 25 1 pixels for sequence II and sequence III 
(resulting in a pixel resolution of 2.0 AU). The images for dif- 
ferent grain sizes are then weighted following the grain size dis- 
tribution and added to compose a final image for each model at 
each wavelength. 

2.3. Comparison to previous simulations 

In the images produced from our simulations (Fig. [T]and[2|), it 
is possible to identify the structures described by (Kuchner "& 
[Holman[ ( [2003| ) using simple geometrical arguments. The gen- 
eral goal and approach o f our simulations is similar to those pre- 
sented by [Stark & Kuchner ( 2008). However, our models cover 
a different region of the parameter space (more massive plan- 
ets, farth er away from the star, la rger dust grains with smaller 
/3). While [ Stark & Kuchner[ ( [2008] ) concentrated on the effects of 
terrestrial planets on exozodiacal dust clouds, the scope of the 
present work is to search for structures induced by giant plan- 
ets in debris disks. Furthermore, we attempt to simulate obser- 
vations of the structures found by our modeling. This requires 
a good sampling of the grain size dependent emissivity of the 
dust and, consequently, of the grain size distribution. The grains 
in St ark & Kuchner (2008) are assumed to be no larger than 
~ 120//m and to be produced by an external source (i.e., a 
planetesimal belt at considerable distance from the modeled re- 
gion), which is a reasonable assumption for exozodiacal dust. In 
the present case, the dust-producing planetesimal belt has to be 
taken into account. Accordingly, larger grains have to be consid- 
ered. Furthermore, these grains cannot be neglected, since we at- 
tempt to simulate observations at (sub-)mm wavelengths, where 
millimeter-sized dust grains significantly contribute to the dust 
emission. Thus, a total of 50 grain size bins is used in the pr esent 
work in contrast to the five bins in St ark & Kuchner [p008| ), and 
we consider grains up to a size of 2 mm. 

Nonetheless, the results from both studies can be compared 
qualitatively. In particular, for models that reproduce similar pa- 
rameter regimes the results are very similar and differences can 
be attributed to the above described differences in the approaches 
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Fig. 1. Simulated images from model sequence I for a face-on orientation of the disk at different wavelengths. The wavelength is 
indicated at the top of each column. The star is at the center of each image. The position of the planet is indicated by the red cross 
in each image. The flux is given in arbitrary units and displayed in a logarithmic stretch from zero to peak value. The region marked 
by a red circle in some images has been attenuated by a factor / to properly display the dynamical range of the images. The value 
of / is then indicated in the lower-left corner of each image. If no value of / is given, no attenuation has been applied. Note that a 
disk without a planet would show no azimuthal structures, but a smooth, featureless radial surface brightness distribution (besides 
inner and outer disk radius). 



and to the diff'erences in the nature of the images themselves. The 
regime of high plane tary mass, far from the s tar in the parameter 
space considered by Stark & Kuchner|2008| can be compared to 
the regime of low planetary mass, close to the star in the param- 
eter space considered in the present work, while images at short 
wavelengths - tracing the small particles in the present work - 
lead to comparable structures. This can be se en from compari- 
son of Fig. 6 (lower left panel) or Fig. 10 of |Stark & Kuchner| 
( |2008j ) to Fig. [T] (model lb at thermal re-emission wavelengths, 
model Id at 24 //m and at 70yum) and Fig. |2] (models Ila and lib 
at wavelengths < 70 jim) of the present work. 

|Wyatt| ([2006) modeled planet-disk interaction in debris 
disks using dynamical simulations including a migrating planet. 
Although our approach is substantially difl'erent, both works at- 
tempt to model structures in debris disks. Therefore, we briefly 
compare the two approaches. In contrast to |Wyatt| ( |2QQ6| ), we 
use no migrating planet to initially trap the planetesimals into 
resonance, but instead assume the dust to be produced through 
collisions in a ringlike, featureless disk. Furthermore, we use a 
much simpler approach to m imic dust crea tion in the planetes- 
imal belt. On the other hand, |Wyatt| ( |2Q06| ) investigated the lib- 
eration of particles from the initial resonances only, while we 
include the possibility that particles are trapped into other reso- 
nances than those in which they are initially (they do not have to 
be in resonance initially at all). Thus, the migrating planet is not 
necessary to produce resonant structures with our approach. 

The difl'erences in the two approaches basica lly result in a n 
inverse situation in the resulting model images. In Wyatt ( |20Q6| , 



very prominent structures are seen at long wavelengths (tracing 
large particles that remain in the resonances which they are in 
initially). At short wavelengths, the disks appear smoother, be- 
cause the particles emitting efficiently at these wavelengths can- 
not be held in the initial resonance. In contrast, our approach 
results in a much smoother appearance of the disk at long wave- 
lengths that trace particles that cannot be trapped into resonance 
efl&ciently, because they are not moving radially through the sys- 
tem. Hence, strong structures at (sub-)mm wavelengths are only 
found in our work if the initial dust distribution is placed close to 
a strong resonance with the planet, or if the dust is close enough 
to the star that even millimeter-sized particles are significantly 
aff'ected by Poynting-Robertson drag. On the other hand, small 
particles are easily trapped into resonances. This is because they 
move radially through the system due to Poynting-Robertson 
drag and, thus, can "find" those resonances. 



3. Results from the dynamical modeling 

From the computed images, one can draw several conclusions 
on the nature and strength of the structures, on the requirements 
for the generation of strong structures, and on the wavelength 
dependence of structures in spatially resolved images. The re- 
sults can be seen in Figs. [T] to [3] and are briefly discussed in the 
following. 
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Fig. 2. Same as Fig.[l] but for model sequences II and III. The flux is given in arbitrary units and displayed in a logarithmic stretch 
from zero to 1/2 of the peak value at wavelengths up to 24 \mv and to twice the peak value at longer wavelengths. Note that a disk 
without a planet would show no azimuthal structures, but a smooth, featureless radial surface brightness distribution (besides inner 
and outer disk radius). 



3A. The face-on case 

For a face-on or nearly face-on oriented disk, the azimuthal and 
radial disk structure can be observed directly. This allows one to 
draw strong conclusions on the disk structure from observations 
directly, without the need of detailed modeling or disentangling 
the real structure from projection efl'ects. On the other hand, it is 
impossible to draw any strong conclusion on the vertical struc- 
ture of the disk. 

Prominent structures at short wavelengths 
Observations at short wavelengths (optical and near-infrared 
scattered light as well as short wavelength thermal re-emission) 
predominantly trace small grains. These grains are subject to ef- 
ficient Poynting-Robertson drag which causes them to move ra- 



dially through the system. Thus, they can easily reach regions 
in the disk, where they can be trapped into resonance by the 
planet. This causes prominent structures in the distribution of 
small grains. The shape and strength of these structures depend 
on the mass of the planet. In general, a more massive planet is 
able to keep a larger region clear of dust than a less massive one. 
On the other hand, less massive planets produce a very promi- 
nent, bar- like structure. 

Prominent structures at small distance from the star 
In a similar way as before, particles at small radii are more ef- 
ficiently aff'ected by Poynting-Robertson drag than particles at 
larger distance from the star. Thus, the same resonant configu- 
ration of planet and disk results in more prominent structures 
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Fig. 3. Same as Figs. [l]and|2l but for an edge-on orientation of the disk. The position of the planet is indicated by the red arrow in 
each image. The flux is given in arbitrary units and displayed in a logarithmic stretch from zero to peak value for model sequence I, 
from zero to 1 / 2 of the peak value at wavelengths up to 24 yum for model sequence II and III, and at longer wavelengths form zero 
to 1 /2 of the peak value for model sequence II and to the peak value for model sequence III. Note that a disk without a planet would 
show a smooth, featureless radial surface brightness distribution (besides the possible signpost of inner and outer disk radius). 



in the distribution of larger grains (at long thermal re-emission 
wavelengths) if placed closer to the star. 

Prominent 1:1 resonance at long wavelength thermal re- 
emission 

Large particles are mostly traced by observations at long wave- 
lengths, because they emit more efficiently at these wavelengths 
and, less important, are cooler than smaller grains. These grains 
are not moving significantly in radial direction through the sys- 
tem because of Poynting-Robertson drag (as long as they are 
not too close to the star, see previous point). Only their eccen- 
tricity causes a small, periodic change of their radial distance 
from the star in addition to the dynamical perturbation by the 
planet. Hence, they can only be trapped into resonance at the 
radial position they are initially placed. This results in promi- 
nent structures caused by 1:1 resonance with the planet, while 
other resonances can be neglected in most cases. The orbital ve- 



locity of these large grains at the same distance from the star 
like the planet is nearly the same as that of the planet due to the 
negligible eff'ect of radiation pressure. This renders the resulting 
structures more stable and, consequently, even more prominent. 
In the case of a narrow, ringlike shape of the disk which is ini- 
tially placed close to the planet, this results in a very prominent 
horseshoe structure (at long wavelength thermal re-emission). If 
the planet is close to the inner edge or within a broad disk, this 
results in a horseshoe-like structure in addition to an otherwise 
nearly featureless disk. If the planet is placed far away from the 
disk, so that it cannot influence the disk by 1 : 1 resonance, the 
eflfect of the planet is very small. In any case, the prominence of 
the structures is increasing with increasing mass of the planet. 

Prominent gaps 

If the planet is placed within a broad disk, it opens a ringlike 
gap in the disk. This gap is very prominent at long wavelengths. 
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while at short wavelengths strong resonant structures dominate 

the appearance of the disk. At long wavelengths, a horseshoe ^ 

structure of large particles trapped into 1 : 1 resonance is visible observing wavelength, 
(see above). The width of the gap and the prominence of the 
horseshoe structure are increasing with increasing mass of the 
planet. 



position of the gaps in the disk. The strength of these wavy struc- 
tures depends on the exact configuration of the system and on the 



3.3. Age dependence of the phenomena 



Faint inner disks 

In all our simulations, we always find a more or less significant 
amount of predominantly small grains which the planet is unable 
to prevent from moving toward the star. These particles move 
inward till they reach their sublimation radius. This inner disk, 
consisting of small, warm grains, is particularly bright at short 
wavelengths (thermal re-emission and scattered light). 

Accumulation of dust close to the planet 
There is a small amount of dust accumulated close to the planet. 
These particles are directly captured by the planet. This re- 
sults in an increased brightness of the planet (and its dust en- 
velope), which increases the chances for direct detection. This 
scenario has als o been suggested to explain the high brightness 
of Fomalhaut b ( |Kalas et al.|2QQ8] ). 

3.2. The edge-on case 

If a disk is seen edge-on, one is faced with a number of chal- 
lenges when deriving its radial, azimuthal, and vertical structure. 
Since debris disks are optically thin, one integrates all flux on the 
line of sight. This results in strong degeneracies between radial 
and azimuthal structures. Furthermore, vertical structures can- 
not be assigned to a particular radial and azimuthal position in 
the disk, but only to a distance from the star projected onto the 
sky plane. However, most of the known, spatially resolved de- 
bris disks are seen close to edge-on. This is most likely because 
of the higher surface brightness due to more emitting material 
on the line of sight which, results in an observational bias. For 
that reason, it is particularly interesting to search for prominent 
structures in edge-on seen debris disks due to planet-disk inter- 
action, which can be unambiguously identified. 

Particles scattered out of the disk midplane 
Particularly for massive planets close to the radial position of the 
dust production, there is a halo of particles scattered out of the 
disk midplane by the interaction with the planet. Therefore, the 
halo consists mostly of small and intermediate- sized grains mov- 
ing efficiently radially through the disk. At short to intermediate 
thermal re-emission wavelengths, resonant structures are visible. 
These structures are particularly prominent at these wavelengths, 
because they mostly consist of intermediate- sized grains. These 
grains are sufficiently affected by Poynting-Robertson drag and 
radiation pressure to interact with the planet and, thus, to be 
scattered out of the disk midplane. On the other hand, Poynting- 
Robertson drag is not strong enough to prevent these grains from 
being trapped into resonance. 

Asymmetries in the disk radial brightness profile 
There is a brightness asymmetry between the two ansae of the 
disk. The disk ansa containing the planet is slightly brighter than 
the opposite ansa. This is due to the accumulation of dust in the 
resonances which results in clumps of higher density (see face- 
on case). 

Multiple peaks and dips in the disk radial brightness profile 
A multi-ring structure as seen in particular in the long wave- 
length results from sequence III results in a wavy radial bright- 
ness profile along the disk midplane in the edge-on case, while 
peaks are seen at the position of rings, and dips are seen at the 



The dominating effect in our model that causes the dust grains to 
change their orbits beside the gravitational interaction with the 
planet is Poynting-Robertson drag. Therefore, one has to com- 
pare the Poynting-Robertson time scale of different grains with 
the age of the systems. This time scale (the time it takes for a par- 
ticle to spiral onto the star starting at a circular orbit at distance 
Rq) can be approximated following (Gustaf son| ( | 1 994j ) : 



400 /M^ 



VAU/ 



yr. 



(2) 



It is short for small grains (large fi) and increases with grain 
size (decreasing JS). It also depends on the radial distance from 
the star. Accordingly, systems with the dust placed closer to the 
star evolve faster. Furthermore, structures in small grains evolve 
faster than large grains. The older a system is, the more time 
there is for larger grains to move significantly in the radial direc- 
tion, and thereby to be trapped into different resonances. With in- 
creasing age of the system, larger grains will have the time to be 
trapped in prominent resonant structures instead of in the single 
1:1 resonant horseshoe structure (in the face-on case). Because 
these larger grains are traced at longer wavelengths, the tran- 
sition between the two types of structures will occur at longer 
wavelengths for older systems. 

It is important to note that no collisions between the dust 
particles are considered in our dynamical simulations. Structures 
that appear only after a long time or that are particularly traced 
by particles that are very abundant may be destroyed by chaotic 
events like collisions. Furthermore, particles may be destroyed 
through collisions before they are able to form these structures. 
The strength of this effect is decreasing with decreasing dust 
mass. Thus, such structures are a result of our modeling, but their 
significance might be over estimated as compared to the strength 
of these structures in real systems (particularly in systems with 
very massive disks). 

3.4. Potential for observations and modeling 

From the above discussions, we found that planet-disk interac- 
tion in debris disks produces structures in the disk that allow one 
to constrain the parameters of a planet-disk system. However, 
data at one of the wavelength regimes considered are expected 
to provide only weak constraints on the actual configuration of 
the system (e.g., planetary mass, major axis, eccentricity, radial 
distribution of the parent bodies), especially when taking into 
account observational effects such as noise, resolution effects, 
and stellar point spread function (PSF) subtraction uncertainties 
(see Sects. |4] to |6] for simulated observations and discussion). In 
the model images themselves, this is particularly obvious, when 
comparing the images at mid-infrared wavelengths (and shorter) 
in face-on orientation produced from our runs II and III (e.g., at 
24yum) and the corresponding (sub-)mm images (Fig. [2]). While 
the mid-infrared data allow one to distinguish particularly well 
between planetary mass and eccentricity (e.g., between run Ila 
and lie), the (sub-)mm data allow one to constrain in particular 
the major axis of the planet and the position of the parent bodies 
(e.g., compare images from runs Ila and Ilia). Thus, combined 
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Table 2. Longest baseline ^max, shortest baseline ^min, and resolution characteristics of the array configurations used for the simu- 
lations. 



Band 10 Band 7 Band 6 Band 5 
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-^min 


-^max 
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^min 


^max 


number 


[m] 


[m] 
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n 
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1.22 
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0.99 


17.2 
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15 
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22.0 
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0.36 


13.1 


0.48 


17.2 


0.61 


22.0 


09 


15 


703 


0.10 


4.8 


0.28 


13.1 


0.37 


17.2 


0.47 


22.0 


11 


15 
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0.07 


4.8 


0.19 


13.1 


0.25 


17.2 


0.32 


22.0 


13 


24 


1440 


0.05 


3.0 


0.14 


8.2 


0.18 


10.7 


0.23 


13.8 


15 


24 


1811 


0.040 


3.0 


0.108 


8.2 


0.140 


10.7 


0.180 


13.8 


17 


49 


2297 


0.031 


1.5 


0.085 


4.0 


0.112 


5.3 


0.144 


6.7 


19 


49 


3105 


0.023 


1.5 


0.063 


4.0 


0.083 


5.3 


0.106 


6.7 


21 


79 


6068 


0.012 


0.9 


0.032 


2.5 


0.042 


3.3 


0.054 


4.2 


23 


79 


11457 


0.006 


0.9 


0.017 


2.5 


0.022 


3.3 


0.029 


4.2 


25,27 


79 


14444 


0.005 


0.9 


0.014 


2.5 


0.018 


3.3 


0.023 


4.2 



observations at both wavelength regimes will be most useful for 
a detailed analysis of planet-disk systems. 



4. Evaluating the observability (1): ALMA 

In the following sections, we simulate observations of the mod- 
eled disk images using diff'erent near- future facilities. In ad- 
dition, we discuss the detectability of the structures with the 
Hubble Space Telescope in order to explain why the prominent 
scattered light structures seen in our simulations have not been 
detected yet. In this section, we start with ALMA. 



4.1. General sensitivity study for debris disks 

ALMA will allow interferometric observations of debris disks in 
the (sub-)mm regime with unprecedented sensitivity and spatial 
resolution. Consequently, it is the most promising instrument to 
study structures in debris disks in this wavelength regime in the 
near future. ALMA is a very complex instrument, providing the 
user with multiple options that allow one to optimize the instru- 
ment set-up for the planned observations. A detailed study of the 
performance of each configuration is necessary to optimize the 
observational outcome and eflftciency. This study in the context 
of spatially resolved observations of debris disk is described in 
the following. 



4.1.1. Model description 

All simulations are carried out for a debris disk around a solar- 
type star (realized as a black body radiator with R^, = LOR©, 
= 1.0 Lq, Teff = 5778 K). A simple, analytical disk model 
is used to keep the results as general as possible. It consists 
of a circular ring (rout = l-lnn) of dust with constant sur- 
face density (similar to the width of a ring of dust particles on 
orbits with identical major axis and eccentricity of 0.1, if the 
directions of the major axes of the orbits are distributed ran- 
domly within the plane of the disk) and an opening angle of 10°. 
Astronomical silicate with a bulk density of 2.7g/cm^ ( Draine 
& Lee||1984t [Weingartner & Draine"2001') is employed for the 
chemical composition of the dust. The grain size distribution fol- 
lows a power-law distribution (exponent -3.5) with lower cut-off' 
size amin = 0.45 yum (i.e., the blow-out size of the system) and 



1.2 
1.1 



a 

i 0.9 



b 

-t-J 

b 
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Fig. 4. Contribution of phase noise to the total noise in our 
ALMA observations (J. P. Ruge, personal communication). The 
noise levels are computed as the standard deviation of all pixels 
in simulated observations of empty images (all pixels set to 0.0). 
The quantity cr^h represents the standard deviation in simulations 
with only thermal noise considered, while (Xth+ph is derived from 
simulations considering thermal and phase noise. Ratios higher 
than 1.0 may occur because of uncertainties in noise estimates 
and because of diff'erent noise contributions that in part cancel 
out each other. The curves for Band 5 and Band 6 are very close 
to 1 .0 and cannot be distinguished in the plot. In any case, the 
contribution of phase noise is lower than 5%. 



upper cut-off* size (2niax = 2 mm. Only observing wavelengths 
up to 1.6 mm are considered. Therefore, the emission of larger 
grains is negligible for a sufl&ciently steep grain size distribution. 
The inner disk radius is set to 5 AU, 50 AU, and 100 AU in dif- 
ferent runs. This co vers a variety of known debris disks, e.g., ou r 
solar system, 6Eri ( Greaves et al.||1998[ Backman et al.||2009|), 
HP 10 7146 dArdila et al.||2004t jCorder eFar r2009 ; Er tel et al. 
'201T), Fomalhaut (Ka las et al.| |2005), and AUMic ( Augereau 
& Beust 2006 ). The dust mass of the disk is set to 10"^ M© for 
simulating the model images and scaled to diff'erent masses later 
(assuming an optically thin debris disk) to evaluate the sensitiv- 
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ity. These masses can be converted into fractional luminosities 
following the equation 



4.5 X 10^ 



AU 



dust 



Mr 



(3) 



i.e., the dust luminosity for a given dust mass and stellar lumi- 
nosity scales with r"^ and is proportional to the dust mass, where 
4.5 X 10^ is the factor for r = 1 AU given the above model of the 
system (stellar properties and disk properties). This is because 
the total luminosity of a given dust species only depends on the 
stellar spectrum, the total surface of the dust (oc dust mass), and 
its absorption efficiency, assuming thermal equilibrium (no in- 
formation about the wavelength range is carried in which the 
majority of this luminosity is emitted, i.e., about the dust temper- 
ature). The resulting factors are 1 .8 x 10^, 1 .8 x 10^, and 4.5 x 10^ 
for the disk models with = 5 AU, 50 AU, and 100 AU in the 
present study. 

Simulations of images from these disk models are performed 
for face-on and edge-on orientation of the disk. The pixel reso- 
lution of the images in AU can then be computed as 2 rout/ 101 
(the total extent of the disk is chosen to be 101 pixels in one di- 
rection) and in arcseconds by dividing the result by the distance 
of the system in pc . Diffe rent distances are assumed to simulate 
observations (Sect. 4.1.2 ). If the total extent of the images as de- 
scribed above is smaller than five times the resolution element 
(FWHM) of the observations (in particular for the models with 
Tin = 5 AU at large distance and small array extent), the back- 
ground fluctuation cannot be quantified in a reliable way. Thus, 
the resolution of the images is lowered in these cases and the 
empty region around the disk image is increased. This is done 
dynamically, so that the pixel resolution is 1/20 of the FWHM 
and the total extent of the images is 200 pixel (10 x FWHM of 
each observation). 

4.1 .2. Simulation of observations 

The CASA ALMA simulatoij^ (procedure simdata) is used to 
simulate observations of the model images. The simulations are 
conducted using the configurations of the ALMA array in full 
operations provided by the CASA simulator (usingonly every 
second configuration from configurations 01 to 21]^ These ar- 
ray configurations contain 50 12-m antennas each. The base- 
line ranges of the array configurations are listed in Table [2] 
Bands 10, 7, 6, and 5 (central wavelengths of 350yum, 950 um, 
1250yum, and 1600yumQare used. The exact central wavelength 
of each bandpass is chosen to avoid strong atmospheric absorp- 
tion bands. The band width used is 7.5 GHz for all simulations 
following the ALMA Cycle Technical handbook^. The ob- 
jects are placed at an optimal position in the sky reaching the 
zenith during observations (RA = IS'^, DEC = -23°). Total 
observing time is 8 h, while the single integration time is 60 s. 
Thermal noise is added to the simulated visibilities. Good, but 
realistic weather conditions are assumed (precipitable water va- 
por = 0.65 mm, ground temperature = 269 Kjj). From a com- 
parison with thermal noise, we find phase noise to be negligi- 
ble (Fig. |4]) and, therefore, neglect this contribution to the total 



^ http://casa.nrao.edu/ 
^ For a detailed explanation 



of all parameters used 
in this procedure, see the corresponding manuals at 

http://casaguides.nrao.e du/index.php ?title=Simulating_ObserYa- | 

tions_in_CASA 

http://almascience.eso.org/document-and-tools 
5 www.eso.org/gen-fac/pubs/astclim/ 



noise. For image reconstruction, natural weighting of the visibil- 
ities is applied. Deconvolution is performed with 500 iterations 
and a threshold of 0.01 mJy/beam. Single pointing observations 
are simulated rather than mosaicking, because most of the simu- 
lated images fit into one field of view of a single pointing obser- 
vation. Mosaicking would require significantly more observing 
time, increasing the total time on target to an unrealistic amount 
considering the high pressure on ALMA expected. 

For the simulated observations of the model images, the 
disks are placed at 1 1 diff'erent distances distributed logarithmi- 
cally from 10 pc to 100 pc. The dust mass (oc total disk flux as- 
suming an optically thin disk) is scaled to difl'erent values. The 
S/N is estimated as the ratio between peak flux in the simulated 
observations of a target and the rms of the background fluctu- 
ation. A linear fit is performed on the distribution of S/N over 
disk mass considering only data with S/N > 20. From this fit, 
the disk mass needed to reach an S/N of 10 is computed. This is 
considered to be a robust detection of the disk. 



4.1.3. Results 

To properly interpret the results, one has to consider the follow- 
ing specifics of an interferometric observation: 

- The field of view of an interferometric observation using one 
single pointing is limited by the FWHM of an observation 
with one antenna (the primary beam). For the 12-m antennas 
of ALMA, this is /I/ 12 m. 

- The spatial resolution ^min of an observation (the FWHM of 
the synthesized beam) is determined by the maximum base- 
line ^max following ^^in = A/B 

max- 

- The largest angular scale ^max of a target that will not be 
filtered out (resolved out) by the interferometric observa- 
tions is determined by the minimum baseline ^min following 
^niax ^ 0.6A/B^^ 

The derived values of ^min and ^max as well as the resolution 
characteristics ^min and ^^ax for the diff'erent array configurations 
and bands considered are listed in Tabled With the derived val- 
ues one can now explain the features seen in the sensitivity maps 
that are shown in Fig. [5] First, there are a number of configu- 
rations in the explored parameter space that cannot be observed 
with the strategy assumed or where the approach used to esti- 
mate the sensitivity produces erroneous results: 

- In particular very extended, nearby disks, cannot be observed 
at short wavelengths without mosaicking, since the spatial 
extent of the disks is larger than the field of view of the obser- 
vations. The configurations where this is the case are marked 
by a black cross in an empty (white) pixel in Fig.|5] Note that 
parts o f the disk can still be observed with one single point- 
ing (Bole y et al.|2Q1 2). 

- For the highest resolution (most extended arrays), in partic- 
ular nearby, extended disks are resolved out. These regions 
are marked in Fig. [5] by an "x". 

- For high spatial resolution, the FWHM of the synthesized 
beam is similar to, or smaller than, the pixel size of the syn- 
thetic images. This results in heavy residuals from image re- 
construction and does not allow one to properly quantify the 
S/N in the simulated maps. Note that this is a limitation of 
the approach, not a real limitation of ALMA. However, at 
such high resolution, the disks would have to be extremely 



^ ALMA Cycle Proposers Guide, |http://almascience.eso.org/| 
document-and-tools 
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Fig. 5. Maps of the 10 cr sensitivity of ALMA to face-on seen (left) and edge-on seen (right) debris disks as a function of the distance 
and array configuration (spatial resolution) for diff'erent wavelengths of observations. The three rows refer to diff'erent disk extents 
given by rin of 5 AU (top), 50 AU (middle), and 100 AU (bottom). Disk mass can be converted into f ract ional l uminosity following 
Eq.[3] For details about the model images and about the simulation of observations see Sects. 4.1.1 and |4. 1.2 The detectable disk 
mass is color-coded and can be converted into fractional luminosity following Eq. [3] White pixels with symbols mark positions 
where the simulations failed or observations are impossible for various reasons denoted by the symbols as follows (see Sect. 4.1.3 
for more details): cross - The spatial extent of the disk is larger than the field of view of a single pointing, x - the disk is "resolved 
out" (the radial width of the ring is larger than the largest scale detectable), open square - artifacts of the image reconstruction 
caused by deficiencies of the model image dominate the simulated observations. 



massive (bright) to be detected at all (see the discussion of 
the sensitivity below). 

In the regions of the parameter space at which the simula- 
tions are not aff'ected by the above limitations, one can evaluate 
the sensitivity. The following behavior in the parameter space 
can be found: 

- In general, the highest sensitivity to disk mass is reached in 
Band 6 (central wavelength 1250yum). 

- For a given spatial resolution, the sensitivity is only slightly 
decreasing with increasing distance as long as the disk is 
spatially resolved (FWHM < 2 rout). This can be explained 
by the fact that the surface brightness of the disk (e.g., in 
mJy/AU^) is decreasing with the distance squared, but the 
beam covers an increasing area of the disk, depending on the 
scales of the dominating disk structures. 

- For a given distance, the S/N is increasing with decreasing 
spatial resolution as long as the disk is spatially resolved. 
This is because the beam covers a larger area of the disk at 
very similar sensitivity (e.g., in mJy/beam). As long as the 
FWHM is smaller than the dominating structure of the disk 
(e.g., the ring width), the sensitivity is increasing with the 
square of the FWHM. 

- As soon as the disk is spatially unresolved, the sensitivity 
drops with the distance squared for a given spatial resolution. 

- In general, disks seen edge-on have a higher surface bright- 
ness. Thus, spatially resolved observations of an edge-on 
seen disk will result in a higher S/N than the same obser- 
vations of the same disk seen face-on. In addition to that, 
there is no significant diff'erence between imaging of edge- 
on and face-on disks (for the employed disk models). This 
can be explained by the fact that for spatially resolved imag- 



ing only the peak S/N is considered, which depends on the 
scale and brightness of the brightest structures in the disk. 
This is in both cases the radial width of the ring and, for disk 
seen edge-on, the vertical height (which is of the same or- 
der in our model). For spatially unresolved observations, the 
orientation of the disk has no eff'ect at all. 

The above qualitative discussion leads to an important recom- 
mendation for optimal spatially resolved imaging of debris disks 
with ALMA. As long as the structures of interest are bright 
enough, one should always use a spatial resolution similar to the 
scales of interest. Most importantly, one should avoid to over- 
resolve the structures, because sensitivity then decreases with 
the square of the spatial resolution. 

ALMA will not only increase our knowledge about the 
known spatially resolved debris disks owing to the higher sensi- 
tivity and spatial resolution. Because of the nearly constant sen- 
sitivity to dust mass in spatially resolved imaging at distances 
up to ^ lOOpc (for disks with radial extents of few tens of AU), 
ALMA will significantly increase the sample of debris disks spa- 
tially resolved at millimeter wavelengths toward more distant 
systems compared to earlier instruments. 

4.2. Observability of planet-disk interaction 

From our dynamical modeling, we select several runs that re- 
sult in prominent structures at wavelengths > 350yum in face-on 
and/or edge-on orientation. The total flux in the selected images 
is scaled to the total flux of selected, well-studied debris disks 
(Table |3]). The CAS A simulator is used to predict results of ob- 
servations. The optimum array configuration is selected by the 
resulting S/N and spatial resolution in the simulated images. 
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Table 3. Reference systems for our simulated ALMA observations 



Run 


Reference 


^ref 


d 


Host star 


Age 


-^dust 


^disk at /iref 


Mdisk 




Ref 






[fim] 


[pc] 




[Myr] 


[AU] 


[mJy] 


[Mo] 






lb 


6Eri 


24 


3.2 


K2V 


850 


3.0 


330.0 


5.7 X 10" 


ll 


1 


Ic 


HD105 


850 


40.0 


GOV 


30 


45 - 120 


10.7 


1.2 X 10" 


-7 


2 


lie 


HD 107146 


350 


29.5 


G2V 


130 


50 - 250 


319.0 


7.3 X 10- 


-7 


3,4 


Ilia 


HD 107146 


350 


29.5 


G2V 


130 


50 - 250 


319.0 


3.7 X 10- 


-7 


3,4 


Illd 


HD 107146 


350 


29.5 


G2V 


130 


50 - 250 


319.0 


5.5 X lo- 


-7 


3,4 



References. (l)|Backman ^taLH2009l); (2)|Nilsson et ari([20T0l); (3)|Corder et aLH2009l); (4)|Ertel et aL|([20TI]). 
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Fig. 6. Top: Simulated ALMA observations of the model image resulting from run Illd (face-on orientation) for different array 
configurations. Only the results from a representative number of array configurations are shown. The number in the upper-right 
corner of each image denotes the array number. The corresponding beam is displayed in the lower-right corner of each image. 
Bottom: Simulated ALMA observations of all selected model images. The optimum array configurations have been selected. The 
run and the array number are displayed in the upper-right corner, the corresponding beam is displayed in the lower-right corner of 
each image. All observations were simulated for Band 6 (central wavelength 1250 yum). 



4.2.1 . Selection of runs and preparation of the images 

The runs for which ALMA simulations are carried out are se- 
lected according to the prominence of the structures seen in the 
(sub-)mm: 

- Runs lb, Ic, lie, and Illd are found to represent the prominent 
structures in the synthetic images well for face-on orientation 
of the disk (Fig. [T| and [2]). 

- Run Ilia is selected as an example for structures seen in edge- 
on orientation of the disk (Fig.jSj. 

From these runs, observations are simulated at Band 6 (central 



wavelength 1250 yum, found in Chapt. [4J] to be the most sensitive 
one for the intended observations) as follows: 

- Known debris disks with radial dust distributions similar to 
the simulated ones are selected as reference objects. Systems 
around solar-type stars are considered, because solar-type 
stars are considered for the dynamical simulations as well. 



- A representative photometric measurement (in the 
(sub-)mm, where possible) is used to scale the total 
flux of the simulated images. Therefore, the total flux of the 
disk at the wavelength of these observations is computed. 

- The disk is placed at the distance of the reference debris disk. 
A flux- scaling factor is derived by comparing the flux of the 
modeled disk with the measured flux at the reference wave- 
length. This scaling factor is applied to the image at Band 6 
to properly scale the flux in this image to a realistic value. 

- The stellar contribution of the host star of the reference de- 
bris disk is added to the center of the synthetic image. 

The reference debris disk for each run, the reference wavelength, 
the applied flux at this wavelength, some basic information about 
the reference system, and the resulting dust mass when scaling 
the model disk to the observed flux are shown in TableO It is im- 
portant to note that the model scaled to the observed flux of the 
reference debris disk at one wavelength is not expected to repro- 
duce the whole spectral energy distribution of the disk properly. 
However, the reference debris disks are selected by their sim- 
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ilarity to the model systems in the radial dust distribution and 
Band 6 is clearly in the Ray leigh- Jeans regime of the dust re- 
emission. Accordingly, the flux at this wavelength is considered 
to be at least realistic for a debris disk similar to the reference 
disk. 

The reference debris disks used for the scaling of the model 
are motivated as follows: 

- 6Eri for run lb: The inner debris ring (~ 3 AU from the star; 
Backman et al.||2009| ) is used. For this ring, the 24 jim flux 



is used as reference value to scale the flux in all model im- 
ages, since the other disk components are expected to have 
significant contribution to the flux at longer wavelengths. 

- HD 105 for run Ic: HD 105 is a massive debris disk around 
a solar-type star. Most of the dust is expected to be con- 
centrated at a distance of 40 AU to 50 AU from the star 
( [Hillenbrand et al.|2008l ). 

- HD 107146 for runs lie. Ilia, and Illd: The work on the 
HD 107146 disk (Ert el et al.||2QTTl ) was the motivation for 
these runs (Sect. |2.l| ). 

4.2.2. Parameters for the CASA simulations 

The simulations are carried out using all configurations of the 
ALMA array in full operations provided by the CASA simula- 
tor. Observing conditions and parameters for the simulations are 
identical to the simulations in Sect. |4.1| Band 6 (central wave- 
length 1250yum) is used. Single pointing observations are simu- 
lated rather than mosaicking, since all simulated images fit into 
one field of view of a single pointing at Band 6. 



image. This can be attributed to the limitations of ALMA to 
observe very extended structures and of the simulations as de- 
scribed in Sect. |4.1.3[ The signal has to be evaluated as the flux 
above this homogeneous background. 

For run Illd, value is placed on the multi-ring structure, 
while the bar-like structure is ignored. An example for an array 
configuration that results in a significant detection of the bar-like 
structure can be seen in Fig. [6] (top). A peak S/N of 10 (array 11, 
FWHM = a.'38 X a.' 47) is found to result in significant detection 
of the three disk components (outer and inner ring as well as in- 
nermost accumulation of dust), since they have a very similar 
surface brightness. The bar-like structure is visible at low S/N 
in these simulated observations, as expected from the previous 
discussion about the connection between spatial resolution and 
sensitivity to surface brightness. 

For the disk from run Ilia seen edge-on, the structure of in- 
terest is the gap in the disk seen as a dip in the radial bright- 
ness distribution and the outer disk seen as a secondary peak 
of the surface brightness beyond this gap. An S/N of 15 (ar- 
ray 16, FWHM = Of! 16 x a'22 - similar to the vertical extent of 
the disk) in the secondary peak of the surface brightness results 
in a clear detection of both features and allows one to marginally 
resolve the vertical extent of the disk. 



These examples illustrate the findings in Sect. 4.1.3 for struc- 
tured disks. ALMA is not able to reach both high sensitivity 
to surface brightness and high spatial resolution simultaneously. 
The highest sensitivity at reasonable spatial resolution is reached 
if the resolution of the observations is similar to the scale of 
the structures expected, as long as these structures are bright 
enough. 



4.2.3. Results 

We find from the simulations of observations at diff'erent spa- 
tial resolution (diff'erent maximum baseline) that a trade-off' be- 
tween sensitivity and spatial resolution is necessary (see also 
Sect. 4.1.3 ). This is particularly well illustrated in the example 
of run Illd (Fig.[6j top). The bar-like structure within the (outer) 
gap in the disk is only visible at the lowest resolutions, while one 
needs a higher resolution to clearly resolve the inner gap and to 
clearly separate the innermost peak from the inner ring of dust at 
V to from the star. For the five model disks considered for the 
ALMA simulations, the best array configurations based on high 
S/N and high resolution are selected. The model image resulting 
from run lb is found to be too faint to be detected by any of the 
simulated observations. For the other four model images, the re- 
sults are shown in Fig. [6] (bottom). The selection of the optimum 
array configuration for each case is described in the following. 
The results are consistent with the suggestions for optimal spa- 



tially resolved imaging of debris disks derived in Sect. 4.1.3 



For run Ic, the structures are very smooth and no substructure 
is seen. Thus, one can use the peak S/N as a good tracer for the 
significance of the structures detected. Based on an additional, 
visible inspection of the images, a peak S/N of 22 (array 19, 
FWHM = a.'ll X a.' 15) is found to give the best results, while 
the next larger array would result in a peak S/N of 15, but large 
parts of the disk would be detected at an S/N < 10. 

For run lie, the structure of interest would be the faint 
horseshoe structure at the inner edge of the disk. This struc- 
ture can in any case only be observed at an S/N of up to 4. 
This is reached when the size of the beam is comparable to the 
size of this structure. The optimum is reached using array 08 
(FWHM = a' 62 X a.'78). It is important to note that there is a 
negative background that is nearly homogeneous over the whole 



5. Evaluating the observability (2): The Hubble 
Space Telescope 

From the dynamical simulations, prominent structures in de- 
bris disks are found, in particular in scattered light and short- 
wavelength thermal re-emission. In contrast, no debris disk seen 
face-on is known to exhibit such structures in scattered light ob- 
servations. Therefore, we brieffy examine the observability of 
such structures in face-on seen disks with present instruments in 
scattered light, in particular with the HST with which most of the 
scattered-light detected debris disks have been discovered. Our 
example here is the HD 107146 debris disk. It is one of the few 
debris disks seen (close to) face-on detected in scattered light 
dArdila et al.|2QQ4l [Ertel et al.|2011| ). 

The HST/ASC image in t he F606W (A^ = 0.6 fim) filter pre- 
sented by |Ardila et al. (^2004 ) is used to evaluate the observabil- 
ity of simulated str uctures in these d ata. Therefore, the model- 
subtracted image of Ertel et al. ( 201l| is employed, which repre- 
sents the pure PSF subtraction residuals (within the capabilities 
of the modeling). To these data, simulated images from the dy- 
namical modeling (inclination of 25° from face-on applied from 
,Ardila et al.,2004) are added after scaling them to the scattered- 
light flux of the fitted disk at the same wavelength. This results in 
simulated observations of the model images using the same strat- 
egy as for the original observations. The results for the model 
images from runs lib and lie are shown in Fig. [7] 

While the disk is clearly detected, the structures of the mod- 
eled images are corrupted by the strong PSF subtraction residu- 
als. Thus, one would not be able to draw strong conclusions on 
the azimuthal structure of the disk from these data. Furthermore, 
there is only a small number of debris disks seen face-on and 
detected in scattered light. In contrast, the two most prominent 
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Table 4. Reference systems for our simulated JWST observa- 
tions 
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Fig. 7. Simulated HST/ACS coronagraphic observations of the 
simulated images from runs lib and lie at 0.6yum. The disk is 
inclined by 25° along the x-axis, so that the lower half of the 
disk points toward the observer. 



debris disks seen edge-on, /SFic and AUMic, exhibit clumpy 
and warped structures that might be interpreted as structures 



Golimowski et al. 


2006|). 


found by the mode 


mg in 



2QQ6| ). All this suggests that the structures 



ing in this work may be present, but not de- 
tected in face-on seen debris disks, yet. 



6. Evaluating the observability (3): Space-based 
near- to mid-infrared telescopes 

Near- future space-based near- to mid-infrared telescopes (e.g., 
the James Webb Space Telescope, JWST, Gardner et al. 2006, or 
the Space Infrared Telescope for Cosmology and Astrophysics, 
SPICA, |Swinyard et al.|[2Q09| ) are expected to provide high- 
sensitivity, high angular resolution imaging capabilities in the 
optical to near-infrared. In this section, observations with such 
facilities are predicted based on available information about the 
expected capabilities of the JWST. Since the focus of this work 
is on spatially resolved imaging capabilities, the two relevant in- 
struments are the Near InfraRed Camera (NIRCam) and the Mid 
InfraRed Instrument (MIRI). NIRCam provides imaging capa- 
bilities in the 0.6yum - 5.0yum wavelength range. In this range, 
coronagraphy is necessary to block the direct stellar radiation. 
The HST is very successful in coronagraphic imaging of bright 
debris disks and the JWST will exceed these capabilities due to 
its higher sensitivity and spatial resolution. However, the perfor- 
mance of coronagraphic instruments depends on several influ- 
ences such as PSF stability. Consequently, the exact results of 
such observations are hard to predict without detailed knowl- 
edge about the performance of the instrument during science 
operations. Hence, the focus of the present section is placed on 
predictions on the capability of debris disk observations through 
none-coronagraphic imaging as are possible with MIRI. 

MIRI will provide imaging capabilities in the wavelength 
range of 5.0 //m - 29yum. In particular at the long wavelength 
edge of this range, thermal radiation of debris disks is expected 
to have a detectable level. With a spatial resolution of V and 
high sensitivity, MIRI is the first instrument that will allow for 
resolved imaging of a large number of debris disks in this wave- 
length range. 



Run 


Reference 


^disk 






Ref. 






[mJy] 


[mJy] 


[Mo] 




la 


6Eri 


330.0 


1726 


8.5 X 10-^' 


1 


lb 


6Eri 


330.0 


1726 


7.3 X 10-12 


1 


Ila 


HD 107146 


19.0 


40.8 


4.4 X 10"^ 


2 


lie 


HD 107146 


19.0 


40.8 


1.4 X 10-^ 


2 


lld 


HD 107146 


19.0 


40.8 


5.7 X 10-^ 


2 


lllc 


HD 107146 


19.0 


40.8 


7.4 X 10-9 


2 



Notes. The reference wavelength at which the fluxes are given is 24 jum 
for both systems. Also note the information given in Tab. |3] about the 
two reference systems. 

References. (l) |BackmanetaLl ( [20091 ); (2) [Hillenbrand et al.| ( [2008l ). 
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Fig. 8. Simulated PSF of JWST/MIRI at 25.5 jim. The image is 
shown in a logarithmic stretch from to 5% of the peak flux to 
highlight the high order structure. The white contours represent 
10% and 50% of the peak flux. 



6.1. Selection of runs and approach to simulate observations 

Several runs are selected from the dynamical simulations to ex- 
plore the capabilities of MIRI to detect the structures. A sum- 
mary of the considered runs and the references used for each run 
is given in Table]?] Since 6Eri is the only known debris disk this 
close and it is seen close to face-on, only the face-on orienta- 
tion is used for the simulated observations on the example of the 
inner disk in this system. For the simulated observations on the 
example of the HD 107146 debris disk, face-on and edge-on ori- 
entations are used, since a number of disks with difl'erent orien- 
tations are known at dista nces of few tens of AU (e.g., A UMic, 
Kalas et al.f2QQ4| Eri, [Liseau et al.||2QlQ| HD 207129, [Krist 
et al.|2QlQ| ). 

Because the star contributes significantly to the flux at mid- 
infrared wavelengths, a stellar PSF subtraction is expected to be 
necessary to clearly reveal any resolved disk structure. To eval- 
uate the observability of structures in debris disks, a two-step 
approach is used: 



(1) Observations of the disk only are simulated by convolving 
the synthetic disk images with the telescope PSF. From the 
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Fig. 9. JWST/MIRI images at 25.5 yum. Leftmost column: Images produced from the runs la and lb. The 6Eri inner debris disk is 
used as reference. Columns two to five: Images produced from the runs Ila, lie, Ild, and IIIc. The HD 107146 debris disk is used as 
reference. The contours in each image represent 3 cr, 5 cr, 10 cr, 20 cr, etc. (double every step) S/N levels for a total integration time 
of 10 s (6Eri) and 900 s (HD 107146). The images are shown in a logarithmic stretch. 



resulting images, the sensitivity needed is estimated and the 
capability to spatially resolve the structures is evaluated. 
(2) The stellar PSF contribution is estimated and the PSF sub- 
traction accuracy needed is evaluated. 

For simulated observations of the disk only, the images at 
24 ^m resulting from the dynamical simulations are scaled to the 
24 yc/m fluxes of two reference debris disks, e Eri and HD 107 146. 
This approach is expected to result in a realistic mid-infrared 
brightness of the disk images, which is important for the sim- 
ulated observations, but not in a good reproduction of the SED 
of the system over the whole wavelength range. The disks are 
placed at the distance of the reference debris disks. The images 
are convolved with a simulated PSF of MIRI at 25.5 jim (Fig.jS]) 
produced by the software WebbPSF^. The resulting sensitivity of 
the observations is estimated using a sensitivity of 30//Jy (10 cr, 
A = 25.5 yum, integration time =10"^ s^and scaling it with the 
square root of the actual integration time ( G ardner et al.|20 06^. 

In Fig. [8] one can clearly see the hexagonal structure and the 
peaks of the PSF at high order. Since the star contributes signifi- 
cantly to the flux in the mid-infrared, one has to perform accurate 
PSF subtraction to distinguish between real disk structures and 
these high order PSF structures. To evaluate the accuracy of stel- 
lar PSF subtraction necessary, the contribution of the PSF de- 
pending on the radial distance from the star has to be computed. 
Therefore, the PSF used to convolve the disk images is scaled to 
the total flux of the star in each of the two reference systems. The 
maximum of the PSF in one-pixel wide radial bins (pixel scale of 
the simulated disk images) is computed. The result is shown in 
Fig. 10 To obtain a reliable image of the disk, the contribution 
of the stellar PSF should be reduced by PSF subtraction (e.g., 
through observations of a reference star) to a level of 1/10 or 
less of the disk flux at a comparable distance from the star. 



6.2. Results 

The results from the simulated observations of the disk only are 
shown in Fig.|9] In both the 6 Eri and the HD 107146 case, one 

^ http://www.stsci.edu/jwst/science/data-simulation-resources; effec- 
tive Aug. 2011 



Table 5. PSF subtraction accuracy necessary for observations of 
the modeled debris disks. 



Run 


Pdisk 

[arcsec] 


^disk 
r mJy j 
[ beam J 


Ppeak 

[arcsec] 


^PSF 
r mJy j 
[ beam J 


cr 

[%] 


la 


1.0 


21.8 


1.4 


90.0 


2.4 


lb 


1.0 


55.6 


1.4 


90.0 


6.2 


Ila (fo) 


2.0 


1.3 


1.4 


2.1 


6.2 


11a (eo) 


2.0 


1.1 


1.4 


2.1 


5.2 


11c (fo) 


2.6 


0.6 


2.3 


0.2 


30.0 


11c (eo) 


3.0 


1.3 


3.2 


0.2 


65.0 


lid (fo) 


3.0 


0.2 


3.2 


0.2 


10.0 


lid (eo) 


3.0 


0.9 


3.2 


0.2 


45.0 



Notes. The quantity pdisk denotes the (projected) distance from the star 
at which the disk flux is measured (the brightest structures of intent), 
Ppeakdenotes the distance from the star of the closest peak of the stellar 
PSlj^ Fdisk and FpsF are the surface brightness of the disk and the PSF 
structures at the corresponding distances, and cr denotes the accuracy at 
which the stellar PSF has to be subtracted to reduce the PSF structures 
to 1/10 of the disk structure (i.e., 1/10 x Fdisk/^psp)- Values are given 
for a face-on (fo) and edge-on (eo) orientation of the disk. 



can see that the disks are spatially resolved. One can clearly dis- 
tinguish between the different planet-disk configurations in both 
the face-on and the edge-on case. Only the results of run IIIc do 
not show a clearly resolved image of the disk. In contrast, the 
bright, unresolved inner part of the disk dominates the emission, 
which results in a bright image of the PSF overlayed on the outer 
disk structure. 

The accuracy of stellar PSF subtraction necessary for the 
diff'erent runs is given in Table [5] This is the result of a compar- 
ison of the simulated disk observations (Fig|9]) and the results 
plotted in Fig. 10 We find that a PSF subtraction accuracy of 
1 % is suflftcient to detect and spatially resolve all simulated de- 
bris disks considered in this study (significantly lower accuracy 
is needed in some cases). Since this uncertainty dominates the 
total uncertainty in the observations, very accurate photometry 
is required to scale the reference PSF to the correct level. This 
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Fig. 10. Radial distribution of the PSF structures at 25.5 //m as 
the maximum flux at a given radial distance. The total flux is 
scaled to the flux of the 6Eri stellar photosphere. To scale it to 
the HD 107146 stellar photosphere, one has to divide it by 42.3. 



represents the minimum requirements to unequivocally detect 
structures in the model debris disks without further considera- 
tions. However, a sophisticated approach might include proper 
rotation of the optics, so that edge-on seen disks are imaged be- 
tween the bright hexagonal wings of the PSF - in regions where 
the PSF structures are less bright - increasing the depth of the 
observations. Furthermore, it might be possible in some cases to 
scale the reference PSF to the high order structures of the PSF in 
regions of the science observations where no significant signal 
from the disk is expected. 

7. Conclusions 

We demonstrated that planet-disk interaction may produce de- 
tectable structures in the dust distribution of debris disks. This 
depends on the configuration of the planet-disk system and 
on the observing wavelength. The detected structures enable 
one to infer and characterize extrasolar planets in a range of 
mass and radial distance from the star unattainable by other 
techniques. In particular, detailed modeling of a combination 
of high-sensitivity, high spatial resolution observations at mid- 
infrared wavelengths and (sub-)mm wavelengths is able to put 
strong constraints on the configuration of the planet-disk sys- 
tem. We demonstrated that HST scattered-light observations are 
in most cases unable to unambiguously detect such structures, in 
particular in debris disks seen face-on. In contrast, both ALMA 
and the JWST will provide the sensitivity and resolution to de- 
tect and spatially resolve the spatial dust distribution in debris 
disks at a level of sensitivity and resolution that allow one to dis- 
tinguish between difl'erent planet-disk configurations. However, 
we also demonstrated limitations of the instruments. ALMA is 
unable to reach both high sensitivity to surface brightness and 
high spatial resolution simultaneously, requiring a sophisticated 
observing strategy to optimize the outcome of planned observa- 
tions. For a debris disk with typical shape of the SED, interme- 
diate observing wavelengths around 1 nmi and small to interme- 
diate array extents are ideally suited to reach a high S/N and rea- 
sonable spatial resolution with ALMA. The highest sensitivity 



^ Of the two peaks of the PSF next to pdisk in radial distance, the 
brighter one is considered. 



at reasonable spatial resolution is reached if the resolution of the 
observations is similar to the scale of the structures expected, as 
long as these structures are bright enough. Mid-infrared obser- 
vations of debris disks with the JWST will be able to detect and 
spatially resolve dust in debris disks even at a distance of several 
tens of AU, where the emission from debris disks in this wave- 
length range is expected to be low. For such observations, stellar 
PSF subtraction with an accuracy of a few percent is necessary 
to unequivocally detect structures in the spatial distribution of 
the dust. 

Acknowledgements. We thank J. P. Ruge for valuable discussions and impor- 
tant contributions to the consideration of phase noise. S. Ertel acknowledges 
financial support from DFG under contract WO 857/7-1 and from the French 
National Research Agency (ANR) for financial support through contract ANR- 
2010 BLAN-0505-01 (EXOZODI), and for general support by K. Ertel. 



References 

Aarseth, S. J. 1985, in Multiple time scales, p. 377 - 418, ed. J. U. Brackbill & 

B.I. Cohen, 377-418 
Apai, D., Janson, M., Moro-Martm, A., et al. 2008, ApJ, 672, 1196 
Ardila, D. R., Golimowski, D. A., Krist, J. E., et al. 2004, ApJ, 617, L147 
Augereau, J.-C. & Beust, H. 2006, A&A, 455, 987 
Backman, D., Marengo, M., Stapelfeldt, K., et al. 2009, ApJ, 690, 1522 
Boley, A. C, Payne, M. J., Corder, S., et al. 2012, ApJ, 750, L21 
Burns, J. A., Lamy, R L., & Soter, S. 1979, Icarus, 40, 1 
Corder, S., Carpenter, J. M., Sargent, A. I., et al. 2009, ApJ, 690, L65 
Dohnanyi, J. S. 1969, J. Geophys. Res., 74, 2531 
Draine, B. T. & Lee, H. M. 1984, ApJ, 285, 89 
Ertel, S., Wolf, S., Metchev, S., et al. 2011, A&A, 533, A132-F 
Gardner, J. R, Mather, J. C, Clampin, M., et al. 2006, Space Sci. Rev., 123, 485 
Golimowski, D. A., Ardila, D. R., Krist, J. E., et al. 2006, AJ, 131, 3109 
Greaves, J. S., Holland, W. S., Moriarty-Schieven, G., et al. 1998, ApJ, 506, 

L133 

Gustafson, B. A. S. 1994, Annual Review of Earth and Planetary Sciences, 22, 
553 

Hillenbrand, L. A., Carpenter, J. M., Kim, J. S., et al. 2008, ApJ, 677, 630 
Holmes, E. K., Dermott, S. P., Gustafson, B. A. S., & Grogan, K. 2003, ApJ, 
597, 1211 

Hughes, A. M., Wilner, D. J., Andrews, S. M., et al. 2011, ArXiv e-prints 

Jewitt, D., Luu, J., & Chen, J. 1996, AJ, 112, 1225 

Kalas, R, Graham, J. R., Chiang, E., et al. 2008, Science, 322, 1345 

Kalas, R, Graham, J. R., & Clampin, M. 2005, Nature, 435, 1067 

Kalas, R, Liu, M. C, & Matthews, B. C. 2004, Science, 303, 1990 

Kokubo, E., Yoshinaga, K., & Makino, J. 1998, MNRAS, 297, 1067 

Krist, J. E., Stapelfeldt, K. R., Bryden, G., et al. 2010, AJ, 140, 1051 

Krivov, A. V. 2010, Research in Astronomy and Astrophysics, 10, 383 

Kuchner, M. J. & Holman, M. J. 2003, ApJ, 588, 1110 

Liou, J.-C, Zook, H. A., & Jackson, A. A. 1995, Icarus, 116, 186 

Liseau, R., Eiroa, C, Fedele, D., et al. 2010, A&A, 518, L132-F 

Liu, M. C. 2004, Science, 305, 1442 

Makino, J. 1991a, PASJ, 43, 859 

Makino, J. 1991b, ApJ, 369, 200 

Makino, J. & Aarseth, S. J. 1992, PASJ, 44, 141 

Marois, C, Macintosh, B., Barman, T., et al. 2008, Science, 322, 1348 

Moro-Martm, A. & Malhotra, R. 2002, AJ, 124, 2305 

Moro-Martm, A. & Malhotra, R. 2003, AJ, 125, 2255 

Nilsson, R., Liseau, R., Brandeker, A., et al. 2010, A&A, 518, A40+ 

Roccataghata, V., Henning, T., Wolf, S., et al. 2009, A&A, 497, 409 

Rodmann, J. 2006, PhD thesis. University of Heidelberg, http://www.ub.uni- 

heidelberg . de/ archiv/ 6123 
Stark, C. C. & Kuchner, M. J. 2008, ApJ, 686, 637 
Stark, C. C. & Kuchner, M. J. 2009, ApJ, 707, 543 

Swinyard, B., Nakagawa, T., Merken, P., et al. 2009, Experimental Astronomy, 
23, 193 

Udry, S. & Mayor, M. 2008, in Astronomical Society of the Pacific Conference 
Series, Vol. 398, Extreme Solar Systems, ed. D. Fischer, F. A. Rasio, 
S. E. Thorsett, & A. Wolszczan, 13-+ 

Vitense, C, Krivov, A. V., & Lohne, T. 2010, A&A, 520, A32+ 

Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 

Wolf, S. & Hillenbrand, L. A. 2003, ApJ, 596, 603 

Wolf, S., Moro-Martm, A., & D'Angelo, G. 2007, Planet. Space Sci., 55, 569 
Wolf, S. & Voshchinnikov, N. V. 2004, Computer Physics Communications, 162, 
113 



15 



Steve Ertel et al.: Observing planet-disk interaction in debris disks 



Wyatt, M. C. 2006, ApJ, 639, 1153 

Wyatt, M. C, Dermott, S. R, Telesco, C. M., et al. 1999, ApJ, 527, 918 
Wyatt, S. P. & Whipple, F. L. 1950, ApJ, 1 1 1, 134 



Appendix A: Orbit integration of mass-less test 
particles 

A.1. Equation of motion 

Because no interaction among the dust grains is considered, one 
can treat them as mass-less test particles independently of each 
other. Summarizing the different force terms, the equation of 
motion governing the dynamics of a dust particle at stellocentric 
position r and velocity v can be written as ([Burns et al.]|1979[ 
ILiou et al.|1995l|Moro-Martm & Malhotra|20Q2| ) 



and its first, second, and third time derivatives to compute the 
time step at which the integration for particle / should be ad- 
vanced in time. The (continuous) time step values are then 
'quantized' to the nearest negative powers of 2 to further econo- 
mize the computation (,MakinO)|1991a| ). Then the computation 
proceeds in blocks of particles that are integrated simultane- 
ously. 

We have tested the accuracy of our numerical integration 
against the semi-analytical solutions of | Wyatt & Whipple| ( |195Q| ) 



for initially circular and elliptical orbit of a dust particle. The 
temporal evolution of the orbital elements and the Poynting- 
Robertson decay time agree very well with Wyatt & Whipple 
(1950 ) . The relative energy error of the orbit integration is typi- 
callybelow 10"^ 



(A.1) 
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where G is the gravitational constant, is the stellar mass, P is 
the number of planets (1 in our case), and Mpi^y and rpi^ are the 
mass and position vector of planet j. The quanti ty is the ratio 
of radiation pressure and stellar gravity force ( Burns et al. | 1979| . 
The strength of the corpuscular stellar wind (consisting mostly 
of protons) is prescribed by the parameter ^, defined as the ra- 
tio of stellar win d drag to Poynting-Robertson drag. Following 
Gustafson ( 1994| ), we adopt the solar value of 0.35. 



/\.2. Orbit Integration 

In this section, we only briefly describe the computati onal pro- 
cedure; for a complete description we refer the reader to"Kokubo| 
et al.| ( |1998] ). The orbits of the dust particles and that of the planet 



are integrated using a fourth-order Hermite scheme, a direct inte- 
gration method using a predict-evaluate-correct (PEC) algorithm 
to sol ve the A/^-body problem ( Makino|1991b(|Makino & Aarseth] 



1992| ). For a particle /, the acceleration aoj and its first derivative 
with respect to time doj are used to predict the future position 
ri(t) and velocity V/(0 from the instantaneous values r^j and Vqj 
at time to < t. 

(1) (3) 

The higher-order derivatives \ and «q ^. are obtained di- 
rectly by Hermite interpolation from a and d at times and t 
( [Makino & Aarseth||1992 ). They are used to improve the ac- 
curacy of the force polynomial to refine the predicted position 
and velocity of the particle. For all simulations we iterate the se- 
quence of (1) evaluation of the acceleration and its first deriva- 
tive and (2) correction of the position and velocity with higher 
derivates three times to increase the accuracy of the orbit inte- 
gration. 

For the choice of the individual time step, we use the pre- 



scription proposed by Aarseth (1985). It uses the acceleration 
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